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We discuss the phase structure of quantum chromodynamics (QCD) with two colors and two 
flavors of light quarks. This is motivated by the increasing interest in the QCD phase diagram as 
follows: (1) The QCD critical point search has been under intensive dispute and its location and 
existence suffer from uncertainty of effective U(1)a symmetry restoration. (2) A new phase called 
quarkyonic matter is drawing theoretical and experimental attention but it is not clear whether it 
can coexist with diquark condensation. We point out that two-color QCD is nontrivial enough to 
contain essential ingredients for (1) and (2) both, and most importantly, is a system without the 
sign problem in numerical simulations on the lattice. We adopt the two-flavor Nambu-Jona-Lasinio 
model extended with the two-color Polyakov loop and make quantitative predictions which can be 
tested by lattice simulations. 
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I. INTRODUCTION 



Understanding the phase structure of quantum chro- 
modynamics (QCD) is one of the key issues in current 
high-energy physics. Thorough phenomenological knowl- 
edge of properties of the hadron spectrum as well as 
nuclear matter is now being complemented by increas- 
ingly precise first-principle numerical studies of QCD at 
nonzero temperature. However, the application of lat- 
tice techniques to matter at high baryon chemical poten- 
tial /is remains a major challenge due to the infamous 
sign problem The difficulties encountered in simula- 
tions of QCD triggered interest in similar theories which 
are free of the sign problem. These include simulations 
at imaginary chemical potential QCD at nonzero 

isospin density p^ - [l6j . a QCD-like theory with adjoint 
quarks pL2-i2C|, and two-color QCD f2ll-[33|. The latter, 
two-color QCD, will be the subject of the present article. 

Two-color QCD differs in several aspects from the 
world we live in. The most notable difference perhaps 
is that the colorless baryons are formed from two quarks, 
and hence are bosons. Dense matter is then not realized 
as an interacting Fermi sea of nucleons, but rather as a 
Bose gas of diquarks which undergoes Bose-Einstein con- 
densation (BEC) at sufficiently low temperature. There- 
fore the ground state of cold and dense two-color quark 
matter forms a superfluid [s^l- Another noteworthy fea- 
ture of two-color QCD, stemming from the fact that 
the SU(2) gauge group has only pseudo-real represen- 
tations, is the Pauli-Giirsey symmetry [35j connecting 
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quarks with antiquarks. As a consequence, the spectrum 
of Nambu-Goldstone (NG) bosons of the spontaneously 
broken chiral symmetry contains diquark states in addi- 
tion to the pseudoscalar mesons. The presence of light 
particles carrying baryon number (i.e. baryonic pions) is 
quite peculiar to two-color QCD (and QCD with adjoint 
quarks) and crucial for understanding the phase structure 
at T 7^ and llb by means of the chiral Lagrangian 
approach [itI . [isl . l36l - [3g | . The model-independent chiral 
Lagrangian arguments have been complemented by in- 
vestigations in various models such as the linear sigma 
model [3^, HO], the random matrix theory j4ll - l43l |. the 
strong-coupling expansion [2^ |43 - |47| . The Nambu- 
Jona-Lasinio (NJL) model was first applied to two-color 
QCD in Ref. 

In fact there has been a tight communication between 
a number of effective model studies and the Monte-Carlo 
simulations on the lattice. Recent progress in this direc- 
tion led to first attempts to probe BEC of diquarks and 
the region of moderate baryon density [29|, H^- Thus, 
based on the knowledge achieved in these preceding pa- 
pers, the present paper aims to make a proposal to use 
the two-color QCD model as a controllable test setting to 
clarify the following controversial issues on the real-QCD 
phase diagram. 

— In-medium U(I)a Symmetry Restoration — 

The QCD Lagrangian has global U(I)a symmetry at the 
classical level. The quantum anomaly, however, leads to 
a nonconserving contribution to the axial current, which 
breaks U(1)a symmetry explicitly. From the point of 
view of quantum field theory the anomaly comes from 
highly ultraviolet modes, and thus, the anomaly should 
be insensitive to any infrared scales such as T, ^lb, and 
the quark mass mg. In this sense the U(I)a anomaly is 
never restored at any T nor /Ltg . In the effective model de- 
scription the U(I)a anomaly manifests itself in the form 
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of a U(l)A-breaking interaction which is microscopically 
induced by instantons i49|] . Because instantons are sup- 
pressed at high T or fiR |5Cll.[5]|. the U(l)A-breaking in- 
teraction is anticipated to weaken in a medium, leading 
to "effective restoration" of the U(1)a symmetry [521 - 
[5^1 • In fact, a quantitative estimate of the reduction of 
the U(1)a effect is crucial for locating the QCD critical 
point in the fiB — T phase diagram [57-59]. The nu- 
merical study by the lattice Monte- Carlo simulation is in 
principle possible in two-color QCD. The clear signal for 
U(1)a restoration is degeneracy in the spectra of mesons 
connected by a U(1)a rotation [HI]. That is, in the two- 
flavor case, the U(1)a partners are 

J Scalar-isoscalar (cr) meson, 

I Pseudoscalar-isoscalar (ryo) meson, 

J Scalar-isovector (oq) meson, 

1 Pseudoscalar-isovector (7?) meson. 

The masses of these multiplets become degenerate when 
the chiral symmetry is also restored. 

Since the a meson involves the so-called disconnected 
diagrams, the lattice simulation is too noisy to see the de- 
generacy with Tyo (flavor-singlet 77) which is also in a noisy 
channel. It is, however, feasible to check the degeneracy 
between ao and 7? in the finite-T and finite- /is lattice 
simulation of two-color QCD. In this paper we will give 
a quantitative guide from an effective model study. 

— Superfluidity and Quarkyonic Matter — 

A new state of matter at high density has been recognized 
and it is now referred to as quarkyonic matter [boI . [61| ] . 
Because the definition of quarkyonic matter is clear only 
in the large Nc limit, in which color non-singlet inter- 
actions are subleading, there seem to be some confusing 
arguments like that quarkyonic matter overwhelms color 
superconductivity. This is not quite true especially when 
Nc is finite. Intuitively, quarkyonic matter is character- 
ized by the properties that the thermodynamic quantities 
(pressure, baryon density, etc.) should be almost sat- 
urated by the degenerated Fermi liquid of quarks and 
the collective excitations on top of the Fermi surface 
should be colorless (i.e. color confined). This is actually 
one of the known properties fulfilled by a certain color- 
superconducting phase, that is, the color-flavor locked 
(CFL) phase [62l465| . To form a colorless object in the 
CFL phase, however, we need to treat a meson composed 
of four quarks, which is technically complicated. Instead 
of real QCD, here, we shall make use of the two-color sys- 
tem; this exhibits superfluidity at high density which is 
reminiscent of color superconductivity in QCD, and be- 
sides, we need not treat four quarks, since diquarks (that 
is, baryons in the two-color world) can make a color sin- 
glet. We will demonstrate that the realization of quarky- 
onic matter is not so exclusive to disfavor superfluidity 
and diquark condensation. The goal of our discussions 
in this part is to make convincing the phase diagram as 
drawn in Fig. [1] 

Let us now specify the model that we will use. Physics 
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FIG. 1: Phase diagram of two-color QCD with two light quark 
flavors. The red and blue bands represent the regions in which 
the Polyakov loop and the normalized chiral condensate, re- 
spectively, take a value ranging from 0.4 to 0.6. The red 
dashed Une extending from the bottom to the right shows the 
onset of diquark condensation. The green band surrounding 
the left-bottom corner indicates the region where ub normal- 
ized by the "Stefan-Boltzmann" (see Sec. IIII Bl for detailed 
explanations) value ranges from 0.4 to 0.6. 



of Cooper pairing of quarks (dense fermionic matter in 
general) is well described by NJL-type effective mod- 
els. However, the NJL-type models suffer from a seri- 
ous drawback: the lack of confinement, which may bring 
about artificial model predictions. In order to capture the 
essential features of confinement physics and yet main- 
tain the technical simplicity of the NJL model, it was 
augmented with the Polyakov loop, equipped with a phe- 
nomenological potential designed to regroduce selected 
lattice data at finite T and zero /z_b O HB-fzl]- (See 
also, Ref. [iS] for example, for a related approach of the 
Polyakov loop extended quark-meson model.) This ex- 
tended model is now called the PNJL model. In spite 
of its simplicity, the PNJL model has shown remarkable 
agreement with thermodynamic quantities measured in 
the lattice QCD simulations. So far, among available 
finite-T and finite- /.ts model options at least, the PNJL 
model is one of the best tools to unveil the QCD phase 
structure. In this paper, in a sense, we downgrade the 
PNJL model into the two-color setting, hoping that our 
model predictions would guide the future two-color sim- 
ulations. 

Figure [1] is actually the phase diagram of two-color 
two-fiavor QCD predicted by means of the PNJL model. 
We note that the phase diagram is divided into three 
regimes: One at high T consists of deconfined and chiral 
symmetric particles, which is to be identified as a quark- 
giuon plasma. The other one at low T and low iib is of 
course the hadronic phase. The last one at low T and 
high /iB is commonly referred to as a superfiuid state. 
We will later discuss that this state can be regarded as 
quarkyonic matter - so to speak, quarkyonic superfluid 
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of two-color quark matter. Keeping in mind this phase 
structure, we will look at the pole and screening masses 
of mesons to extract the information on U(1)a symmetry. 
The meson spectrum in the PNJL model was first ana- 
lyzed in Ref. [69]. Along the same line, we will perform 
the calculations and examine the dependence of meson 
spectrum on the U(l)A-breaking interaction strength. 

The paper is organized as follows. In Sec. we in- 
troduce the model Lagrangian and derive some basic an- 
alytic formulas. Before going into the numerical study, 
we describe in Sec. IIII Al in detail the way we fix the pa- 
rameters of our model. Most of the results have been 
obtained numerically and are presented in Sec. IIIII Fi- 
nally, in Sec. IIVI we summarize and conclude. 



II. MODEL SETUP 

Two-color QCD with N f massless quark flavors has 
a global U(2iV/) flavor invariance at the classical level 
owing to the Pauli-Giirsey symmetry [35| . The axial 
anomaly explicitly breaks U(2A^/) to SU(2iV/). In the 
vacuum, the SU(2iV^) symmetry is spontaneously broken 
by the standard chiral condensate down to its Sp(2A^/) 
subgroup. In the Nj = 2 case, which is the subject of 
the present paper, the symmetry-breakin g pa ttern can be 
equivalently cast as S0(6) S0(5) [13, iHl. The spec- 
trum of NG modes therefore consists of a single 5-plet, 
including three pions and two diquarks (a diquark and 
an antidiquark - a baryonic and an antibaryonic pion). 
Within the NJL model, this degeneracy is reflected by 
the equality of couplings in the meson and diquark chan- 
nels |48| . 

Let us begin the NJL analysis with the U(l)A-invariant 
interaction Lagrangian, 

where a and r denote Pauli matrices in color and fla- 
vor spaces, and ip'^ the charge conjugation of the Dirac 

spinor — C?/' with C — 17^7'^. The interaction Ci 
is minimal in the sense that it only involves the scalar 
and pseudoscalar channels with isospin zero and unity. 
While the above £1 is invariant under a U(1)a rotation, 
we further need an interaction which breaks the U(1)a 
symmetry. To that end we consider the analogous inter- 
action as follows: 



aG 



+ |V'^750'2T2V'P - |V'^cr2T2-0|^ • (2) 



The general interaction Lagrangian is thus a sum of these 
two; Cint = Ci + £2- Ata = only Ci remains and the 
interaction preserves U(1)a, whereas at a = 1 the re- 
maining piece C2 breaks U(1)a maximally, being equiva- 
lent to the two- flavor instanton-induced interaction |3J| . 



Indeed a is a usually used U(l)A-violating parameter [TJI 
but we will also use C = 1 — 2q! for notation simplicity. 
In previous works, ^ = (a — 1/2) was used to discuss 
the phase structure [i^, [t^] • 

Performing the Hubbard-Stratonovich transformation 
in all six channels, we arrive at the total Lagrangian as 

C = V(*7^^M^"^o + 7o/-i-a--j75'?-7?-jC75?7-C'?-a)V' 
+ i (^A*^ij5a2T2ij + h.c.) + ^ (^AfpCia2T2i^ + h. 

- ^(^^ + + Cv' + Ca' + |Ap + CI . (3) 

The covariant derivative Z?^ involves coupling of the 
quarks to the background gauge field A4 which trans- 
lates to the Polyakov loop in the end. In addition, 
we note that mo and /i denote the current quark mass 
and the quark chemical potential. Later we will in- 
troduce iiB to denote the baryon chemical potential; 

= 2ijL where 2 comes from the number of colors. 
The collective fields a, tt , 77, a. A, A5 represent in order the 
mesons in the scalar-isoscalar, pseudoscalar-isovector, 
pseudoscalar-isoscalar, and scalar-isovector channels, and 
the scalar and pseudoscalar diquarks. (We hereafter omit 
the subscript "0" out of 770 and ao for simplicity.) 

In the absence of isospin chemical potential the isovec- 
tor modes do not develop a vacuum expectation value. 
Moreover, the Vafa-Witten theorem [t^ guarantees that 
the chiral condensate in the vacuum has positive par- 
ity. In our model approach, this requires a > 0, that 
is, C ^ 1- (At the same time, C ^ 0> i-G-, a < 1/2 is 
needed for our mean-field treatment using the Hubbard- 
Stratonovich transformation [t^-) Consequently, with 
the exception of the U(l)A-conserving limit a = 0, the 
scalar chiral and diquark condensates will always be pre- 
ferred to the pseudoscalar ones. We will therefore take 
into account only the a and A condensates. The mean- 
field thermodynamics of the system is then independent 
of the parameter a, which will only affect the propagation 
of collective modes, to be discussed in Sec. Ill Bl 



A. Thermodynamics 

In the PNJL model, one introduces a constant tempo- 
ral gauge field which couples to the quarks via the covari- 
ant derivative. In the Polyakov gauge this gauge field is 
diagonal in the color space, and for the color SU(2) group 
it has a form, = iAq = t^O^ where is a real "phase". 
The order parameter for deconfincmcnt is then the traced 
Polyakov loop given by 
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(4) 



where /3 is the inverse temperature. In the mean-field 
approximation, the thermodynamic potential is given by 
a sum of the gauge and quark parts. 



quark ■ 



(5) 
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In the following, we will refer to the two quark colors for 
simplicity as the "red" and "green" . Combining the red 
quark and the green antiquark into the Nambu-Gor'kov 
spinor, ^ = {■ipr,T2ipg)'^ , the background gauge field be- 
comes proportional to the unit matrix in this doubled 
space and the quark thermodynamic potential can be ex- 
pressed as 



quark 



AG 



(fik 
(2^ 



Tr log(zw„ -i0- Hk) , (6) 



where the Nambu-Gor'kov Hamiltonian reads 
^a k + 7oM - ^ -7075A 



7075A* ak + jqM + /i 



(7) 



Here AI — mo + cr is the constituent quark mass and the 
a and A now stand for the condensates. The four eigen- 
values of the Hamiltonian are easily found as and 
—E^ corresponding to the gapped quasiparticle disper- 
sion relations, where 



A2 



The total thermodynamic potential thus becomes 



(8) 



n = -hT[2A^'^e~f^°- + log(l - $2)] + 
(Pk 



AG 



i=± 



(9) 



The first term is the gauge part figaugo having two model 
parameters a and h. We assume the simple form moti- 
vated by lattice strong-coupling expansion [s^ . It differs 
from the three-color expression by a simpler logarithmic 
term due to the SU(2) Haar measure, and by the rescaled 
prefactor of the exponential term, which is in general pro- 
portional to N^. Note that, as usual in the PNJL model 
literature, we simulate the effects of gauge dynamics by 
a constant background temporal gauge field. We then 
adopt a phenomenological ansatz for the gauge contribu- 
tion to the mean-field thermodynamic potential, chosen 
to reproduce selected features of the pure gauge theory. 
Therefore, the parameters a, h only enter the thermody- 
namic potential © since there are no dynamical gauge 
degrees of freedom in our model Lagrangian (|3]). 

It is interesting to recall that in the three-color case the 
thermodynamic potential in general cannot be written in 
terms of the Polyakov loop variable $ (and the conju- 
gate $) only, and one has to use two phases analogous 
to our 9 to parameterize it. On the contrary, in the two- 
color PNJL model the thermodynamic potential depends 
just on $ even in the presence of a diquark condensate. 



This is because diquarks are colorless. This consider- 
ably simplifies the discussion and also avoids technical 
ambiguities stemming from generally complex effective 
actions involving the diquark condensate [H, [7l|, [zl, [tI] . 
The values of the condensates in thermodynamic equilib- 
rium are determined by minimizing the thermodynamic 
potential with respect to the variables cr. A, and $. 

In the quark sector, the effect of the Polyakov loop 
as compared to the simple NJL model is to make the 
replacement E + 2riog(l -I- e'^^) E + riog(l + 
2^e~^^ + e^^^^) in the quasiparticle contribution to the 
thermodynamic potential. Similarly, in the gap equations 
as well as collective mode propagators, one generalizes 
1 - 2f{E) = tanh(/3£;/2), where f{E) = l/(e^^ + 1) is 
the Fermi-Dirac distribution, to 



smhilSE) 
cosh{l3E) + $ 
1 - $ 



1 



cosh(/3i;) -I- $ 
$ 

cosh{l3E) + $ 



[1 - 2fiE)] 
[l-2f{2E)]. 



(10) 



(11) 



The latter two forms ((TU)) and (jTT]) of the function (f{x) 
illustrate that for $ — )■ 1 thermal excitations are domi- 
nated by quark modes with baryon number 1/2 (i.e. de- 
confinement), while for $ — > by baryon modes with 
baryon number 1 (i.e. confinement). With these replace- 
ments, one can readily generalize the results of the pure 
NJL model to the Polyakov loop-extended one, as ob- 
served in Ref. [g^ for the three-color case. 

At zero temperature the Polyakov loop expectation 
value is zero for all values of the chemical potential. This 
led to the suggestion that the PNJL model can naturally 
describe quarkyonic matter [gO, HH at high chemical po- 
tential, that is, a phase where chiral symmetry is restored 
but confinement persists [tI, [7§| . A simple glance at 
Eq. ^ shows that when $ = 0, thermal excitations are 
indeed governed by the term e"^^^*", i.e., they corre- 
spond to colorless baryons. Let us take a closer look at 
how $ ~ arises as T — 0. 

We consider the gap equation for $ following from 
Eq. dH), that is. 



5$ 



1 



1 - $2 



24e" 



Pa \ _ 



E 



d^k 



{2t:)3 ^ + cosh{f3El) 



(12) 

At high enough /i and low T the system is in the Bardeen- 
Cooper-Schrieffer (BCS) regime where Cooper pairing of 
quarks occurs close to the Fermi sea [t^ . The right-hand 
side of this equation is then dominated by the particle 
{i = —) part. We can further simplify the calculation 
by using the high-density approximation, in which we 
expand the dispersion relation around the Fermi surface, 
E^ — _|_ ^ A-|-^2/(2A), and replace the measure 
(fik/{2TTY by A/'dC: where U = fikp /{2Tr^) is the density 
of states at the Fermi surface. The integral thus becomes 
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Gaussian near T ~ and we arrive at the asymptotic 
result as 

*Bcs«^e-''^V2^. (13) 

The Polyakov loop itself is hence suppressed exponen- 
tially at low temperature. 

B. Collective modes 



where w„ = w„ — and ujn stands for the fermionic Mat- 
subara frequencies. The constant Ci is equal to 1/(2G) 
for cr,7r, to 1/(4G) for A, A*, to C/{2G) for 77,0, and to 
C/(4G) for A5, Ag. The trace is taken in Dirac, flavor, 
as well as Nambu-Gor'kov space. When calculating the 
partial derivatives of the Hamiltonian, that determine 
the Yukawa couplings of quarks to the bosonic modes, 
we demand that all fields have to be kept as in Eq. 



The propagators of the collective modes are easily ob- 
tained by a second variation of the effective action that 
follows from Eq. ^ after integrating the quarks fields 
out. Denoting a set of collective fields symbolically as 
Xi, the inverse propagator at imaginary (bosonic Mat- 
subara) frequency iw^ is given by 



, r fj'iu 



X Tr 



dU 



(27r)3 

dn 1 



In the diquark condensation phase, the baryon number 
is spontaneously broken. As a result some of the modes 
mix and we have to find their dispersion relations by 
diagonalizing a matrix propagator. This applies to the 
scalar-isoscalar modes, a, A, A*, as well as pseudoscalar- 
isoscalar modes, 77, A5, Ag. The only modes that do not 
mix are 7? and a thanks to conservation of isospin and 
parity. Their inverse propagators are trivial in isospin 



(14) space, namely, 



S^JD- 



After analytic continua- 



tion to real frequencies they acquire the following forms: 



1 

'2G 



y y f- 



OJ - 



q+f q- 



'2G 



i,j=± k,l=± 



^ q+f q-f ■' ^q- 

oj + kE\.p + IE~\ 




(15) 



Here, we have changed the momentum notation from A; to q to reserve k to take a summation over ±. It should be 
noted that in the limit |p| — > 0, only the i = j terms survive and the expressions somewhat simplify. The first of the 
formulas can be used together with the gap equation for A to show that in the diquark condensation phase the pion 
(pole) mass is exactly equal to 2/i = /is at zero temperature. The proof is straightforward, hence we omit the details. 



In the normal phase (where A = 0) all propagators can be evaluated easily. We provide here the list of expressions 
that we later in Sec. IIII Cl use to calculate the masses numerically. To recall that these are the propagators in normal 
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matter, let us write a superscript (n) 



r)(«)-i 



r)(")-i 
n 



D 



(")-i 



D 



1 

2G 



1 

2G 



2G 

2G 
1 

AG 

S_ 

4G 




ijj 



l + ij 




r\> — p ^ ri, ^ 



+ + jefc_P 



(16) 



Propagators of the antidiquarks, A* and Ag, are ob- 
tained from the diquark ones by charge conjugation, i.e. 
changing the sign of the chemical potential. 



III. NUMERICAL RESULTS 

Armed with these analytic expressions we are now 
ready to proceed in the numerical calculations. First we 
will explain our choice of the model parameters. Next we 
will investigate the phase diagram, in which we particu- 
larly look into the nature of the superfluid phase. Finally 
we will come to the discussion of the meson spectrum de- 
pendence on the U(l)A-breaking parameter a. 



A. Parameter fixing 

Our model has five parameters whose values have to 
be chosen appropriately: a and h in the Polyakov loop 
sector, and G, toq, and the sharp three-momentum cutoff 
A in the NJL sector. (The parameter a or C will be 
treated as an undetermined free parameter.) 

Let us first concentrate on the latter. In the three- 
color NJL model one normally fixes the values of G, mo, 
and A from the physical pion mass rriTr, decay constant 
/tt , which are both measured experimentally, and the chi- 
ral condensate, which is calculated on the lattice or es- 
timated from the QCD sum rules 7J]. The authors of 
Ref. [il] chose to set up the two-color NJL model using 



roughly the same input values, while in Ref. [75| the pion 
mass, decay constant and the (three-color) constituent 
quark mass were used, resulting in a rather different pa- 
rameter set. We shall argue here that both these fits may 



be actually overdetermined and one should be very care- 
ful especially when comparing the model outputs with 
lattice data. 

The chiral limit QCD, as a Yang-Mills theory coupled 
to massless quarks in the color fundamental representa- 
tion, has a single mass scale, namely Aqcd- All physical 
quantities should be expressed as an appropriate power of 
Aqcd times a dimensionless number. In the NJL model 
one can therefore choose freely a single input quantity, 
basically just to set the energy unit. Dimensionless com- 
binations of observables are a pure prediction of QCD, 
and shall thus not be tuned to arbitrary values. In par- 
ticular, either the pion decay constant, the chiral con- 
densate, or the constituent quark mass may serve for this 
purpose, but not two of them independently. With the 
fact in mind that the current quark mass is a free param- 
eter in the lattice simulation, we note that the pion mass 
can in principle acquire any value and represents a new 
physical scale in addition to Aqcd- 

Unfortunately, we are not aware of any lattice data 
which would provide us with the input needed to fix the 
NJL model parameters unambiguously. In order to make 
at least an educated guess, we use an argument based 
on the Nc scaling of physical quantities. Using the fact 
that is proportional to ^/W^ and the chiral condensate 
to Nc, we rescale the three-color values by factors 
and 2/3, respectively. Regarding the Polyakov loop sec- 
tor parameters, the constant a is related to the criti- 
cal temperature Tc for deconfinement in the pure gauge 
theory by a = Tc log 24. Since Tc in the first approx- 
imation does not depend on Nc [60|, we use the value 
Tc = 270 McV. The parameter h can in principle be 
adjusted in order to make the chiral and deconfinement 
crossovers happen at about the same temperature [ssj . 
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TABLE I: First two lines: physical quantities used as an input. Last two lines: fitted parameters of the model. 



Tc [MeV] 


[MeV] 


-(V>„^„>i/3 [MeV] 


[MeV] 


[MeV] 


270 


425 


218 


75.4 


140 


a [McV] 


h^/^ [MeV] 


A [MeV] 


G [GeV-2] 


mo [MeV] 


858.1 


210.5 


657 


7.23 


5.4 



o 
b 



0.0 




12 3 

FIG. 2: Chiral and deconfinement crossovers at /^b = 0. The 
thick lines show the PNJL result (black solid: chiral conden- 
sate in units of the vacuum value ao; red dashed: Polyakov 
loop <1>), while the thin lines indicate the pure chiral (NJL) 
and pure gauge transitions. As expected, the deconfinement 
transition in the absence of quarks is of second order. 



Here we use the estimate based on lattice strong-coupling 
expansion, h — {cs/af' , where as — (425 MeV)^ is the 
physical string tension. The input values as well as the 
fitted parameter set are summarized in Table HI 

In order to check that our parameter set is reason- 
able, we plot in Fig. [2] the expectation values of a and 
$ as a function of T at ~ 0. The positions of the 
two crossovers move very close to each other when the 
coupling between the quark and Polyakov loop sectors is 
switched on. We note that this observation of simultane- 
ous crossovers is less clear if we use the unrescaled input 
parameters. 



B. Phase diagram and diquark condensation 

In the vacuum the diquark is degenerate with the pions 
and its mass is therefore m^r. Thus, when /is(= 2^) ex- 
ceeds TOtt , the diquark condenses and one enters the BEC 
phase which forms a superfluid component. At this mo- 
ment the constituent quark mass M is still rather large. 
However, as further increases, M drops. Once ^ > M, 
a Fermi sea of quarks appears. Here we note that the 
baryon number density, ng, acquires a contribution from 
the diquark condensate and becomes nonzero as soon as 
Hb = TOtt or = TO7r/2. Because of the binding energy, 



naturally, is smaller than M , and thus the onset of 

ns 7^ emerges first and then a quark Fermi sea shows 
up with increasing ^. In the presence of the Fermi sea 
the diquark condensation is closer to the BCS pairing of 
quarks sitting near the Fermi surface rather than to BEC 
of bound bosonic molecules. There is no phase transition 
associated with this qualitative change of behavior, so 
one speaks of a BCS-BEC crossover. Even though it is 
not particularly sharp, it can be conveniently defined by 
the condition /i = M [23 . 

We will draw the phase diagram of our model in a 
"conventional" way. In the phase diagram the deconfine- 
ment crossover is conveniently defined by the condition 
$ = 0.5. The deconfinement temperature is then almost 
independent of the chemical potential. In fact, the value 
of the Polyakov loop at all temperatures depends on /is 
very weakly, as can be seen from Fig. |31 where all con- 
densates are plotted as a function of T/ttIt^ and /is /ttLt^ 
where is fixed to be 140 MeV. This behavior of the 
Polyakov loop can be traced back to the fact that the 
two-color diquark is a color singlet, and therefore does 
not break center symmetry to induce nonzero Polyakov 
loop directly. 

Even at nonzero quark mass the diquark condensate 
exhibits a clear second-order phase transition as shown 
in the plot for A/ctq (where ctq is the vacuum value of the 
chiral condensate) in Fig. [3] Unlike the deconfinement 
crossover, we can draw a well-defined phase boundary in 
the phase diagram separating the normal and superfluid 
phases. As we have already mentioned above, we can 
confirm that nonzero A certainly appears at /is = at 
T = 0. 

A compilation of these data leads to the phase diagram 
we present in Fig. 3] The solid line represents a second- 
order phase transition between the normal and superfluid 
phases. In the superfluid region we have added a dash- 
dotted line which indicates /i = M and can be interpreted 
as a BCS-BEC-type crossover. The dashed line is the de- 
conflnement crossover deflned by = 0.5. Figure [4] is a 
basis to understand a more "advocative" way of present- 
ing the phase diagram as shown in Fig. [1] 

Now we are well prepared to discuss the physical mean- 
ing of each phase labeled in Fig. [T] The red band which 
spreads almost straight along the horizontal axis repre- 
sents the deconfinement crossover. Because a crossover 
has a width and does not have a unique definition, it 
should be much more reasonable to express the transi- 
tion region not by a line but a band. The band width 
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A/ao 




FIG. 3: Chiral condensate a, diquark condensate A, which are given as divided by the chiral condensate in the vacuum ao, 
and the Polyakov loop $ as a function of and T in units of the pion mass = 140 MeV. Note the orientation of the axes, 
chosen to obtain a better view of the surfaces! 
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FIG. 4: Conventional presentation of the phase diagram of 
two-color QCD from the PNJL model in the hb — T plane. 
We indicate the chiral-symmetry breaking (xSB), BEC, and 
BCS phases. Solid line: critical onset of diquark condensation 
where A starts to be nonzero. Dashed line: deconfinement 
crossover (defined by $ = 0.5). Dash-dotted line: BCS-BEC 
crossover (defined by ^ = Af). 



in fact tells us how rapid or slow the crossover is. We 
drew the deconfinement band by the condition that $ 
ranges from 0.4 to 0.6. The blue band showing a sud- 
den decrease around hb/itItt ~ 1.5 is the chiral crossover 
defined similarly by the condition that a/ao ranges from 
0.4 to 0.6. We see that the two crossovers of deconfine- 
ment and chiral restoration take place simultaneously at 
zero density, and this coincidence persists until around 
fJ-B /fTT-n — 1.5. The pink dotted line represents the super- 
fluid onset, which is a well-defined phase transition. Fi- 
nally the green band which shows behavior similar to the 
chiral crossover is drawn by the baryon number density 
n B ■ We can compute hb in the PNJL model and normal- 
ize it by the "Stefan-Boltzmann" value. Since the PNJL 
model is a cutoff theory, even the free (non-interacting) 
limit suffers from the cutoff artifact and deviates from 
the standard formula, which is also the case in the lat- 
tice simulation [80| . Therefore we evaluate the baryon 



number density in the Stefan-Boltzmann limit (nB)sB 
using the PNJL model with M — and $ = 1 imposed 
by hand. In this way, we indicate by the green band in 
Fig.[l]the region in which nB/{nB)sB ranges from 0.4 to 
0.6. 

Because the chiral phase transition controls the dy- 
namical quark mass, it is conceivable that the BCS-BEC 
crossover in Fig.|3]is associated with the chiral crossover. 
This is indeed the case; the dash-dotted line in Fig. S] is 
covered by the chiral crossover band in Fig. [T] Hence, 
as labeled in Fig. [T] the right-bottom region is charac- 
terized by small $ (confinement) and small a/ao (chiral 
symmetric) . This is in fact in accord with the identifica- 
tion of quarkyonic matter in Refs. [H, [t^. In this case 
the blue band is interpreted as the quarkyonic transition. 
If we use ub/ (»-_b)sb to define the quarkyonic transition 
according to Refs. [60,|81[, the green band, instead of the 
blue band, plays the role of the quarkyonic boundary. 

We would emphasize here that the former criterion 
makes more physical sense at least in the present case 
of Nc = 2. Actually, our consideration of the BCS-BEC 
crossover provides us with a clear view point on this issue. 
As we have already discussed, in the chemical potential 
window 1 ^ /is/mTT ^ 1-5, finite baryon number density 
grows. The carriers of the baryon number are, however, 
not quarks but baryons (baryonic pious). This is so be- 
cause < Ub < 2M in this region. Therefore, in the 
phase region we called BEC in Fig.|4l the more appropri- 
ate physical interpretation should be "superfiuid nuclear 
matter" rather than a quarkyonic state in which the pres- 
sure is mostly given by Fermi-degenerated quarks. This 
difference in the interpretation makes a contrast to the 
large- A^c arguments .60, 81] , and it still remains a ques- 
tion to which the real world with Nc = 3 is closer, infinite 
Nc or Nc = 2? Answering this question goes beyond our 
current scope. If the quark-hadron continuity scenario 
driven by the CFL state is realistic, we can say that the 
situation at Nc = 2 is more relevant. 

Before closing this section we mention the previous 
studies. In preceding works there was some controversy 
regarding the order of the phase transition from the di- 
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a 

FIG. 5: Dependence of the vacuum meson masses on the 
U(l)A-breakmg parameter a. The black sohd is m-n-, the red 
dashed rua, the blue dotted m,,, and the green dash-dotted 
ma- We see that m,r = and raa — ma in the U(1)a 
symmetric (a = 0) case. 
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FIG. 6: Meson real-part masses at hb = as a function of 
T/m^. The notation for the lines is the same as in Fig.[Sl The 
black solid is m-^, the red dashed m^, the blue dotted m^, and 
the green dash-dotted ma- The thick lines correspond to 
and ma at a = 0.05, while the thin lines to a = 0.1. 



quark condensation phase to the normal phase. In |48[ it 
was concluded that there is a tricritical point at /.ts/m^ 
somewhere in the range 2.2 — 2.4, and for higher fis the 
transition becomes first order. On the other hand, the 
authors of [t^ used the Thouless criterion to calculate 
the critical temperature, which assumes that the transi- 
tion is second order. Our numerical results for the di- 
quark condensate as a function of T/m,r and fJ-B/m-jr 
(see Fig. |3]) suggest that the transition is second order 
everywhere. We have further confirmed the statement 
that there is no first-order phase transition by looking at 
the quartic Ginzburg-Landau coefficient in our numer- 
ical calculation. Detailed computations and arguments 
are given in the Appendix. 



C. Collective mode spectrum 

We are going to investigate the mass spectrum of col- 
lective modes as a function of T and fi b for different val- 
ues of the U(l)A-breaking parameter a. A useful starting 
point therefore is the a-dependence of the masses in the 
vacuum. This is shown in Fig. [SJ It is worth noting that 
we only display the meson spectrum; thanks to the un- 
broken SO (5) symmetry at = 0, the scalar diquark 
A is degenerate with tt mesons and the pseudoscalar di- 
quark A5 is degenerate with the ao mesons. We should 
perhaps emphasize that whenever we speak of a (pole) 
mass, we have in mind the zero of the real part of the in- 
verse propagator [s^ . We will therefore sometimes refer 
to it as the real-part mass. This coincides with the posi- 
tion of the pole if the zero appears below the threshold 
for decay into quark pairs. Another prescription for the 
mass would be to take the real part of the complex pole 
of the meson propagator, yielding a somewhat different 
result, or to compute the spectral function whose peak 
position and broadness indicate the physical mass and 



decay width. 

In the U(l)A-symmetric limit (a — 0), in our calcula- 
tions, the masses of rjo and ao are equal to those of tt and 
(7, respectively. The degeneracy of rjo and tt is generally 
exact only when the quark mass is strictly zero, so that 
they are both massless NG bosons. In the (P)NJL model 
in the mean-field approximation, this exact degeneracy 
holds regardless of finite quark masses, which is an arti- 
fact of the approximation. In fact, 770 and tt belong to dif- 
ferent irreducible representations of the unbroken SO (5); 
770 is a singlet and tt (with A) form a quintet. The de- 
generacy is only approximate for small quark mass once 
higher-order meson loops are taken into account. 

Off the limit of a — 0, both rriri and rua increase 
steeply. The typical value of a in the three-color NJL 
model with two light quark flavors, that one can obtain 
either by fitting the physical 77' mass or by a reduction 
of the three-flavor model, is in the range 0.1 — 0.2 [t^ . 
In order to be able to investigate the convergence to the 
U(l)A-symmetric limit, we consider in the following two 
particular values, a — 0.05 and a = 0.1. That is, if we 
assume that the vacuum value is a ~ 0.1, we consider 
two examples of no reduction at all and 50% reduction 
of U(1)a effects. It is quite unlikely that a goes to zero 
in the T and /is range of our interest. 

In Fig. [5] the meson real-part masses are plotted as a 
function of T/ttt,^ while keeping fiB = 0- The reason 
why some curves exhibit a cusp structure is that the cor- 
responding modes cross the quark-antiquark (or quark- 
quark in the case of diquarks) threshold as T increases 
(and M decreases accordingly). Of course, such a decay 
into quarks is unphysical and is just an artifact, following 
from the lack of confinement in the NJL model. Even in 
the PNJL model these decay processes are not sufficiently 
suppressed [i^, since the coupling to the Polyakov loop 
only imitates confinement in a statistical sense. 

The real-part masses are difficult to measure on the lat- 
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FIG. 7: Meson screening masses at [Ib = as a function of 
T /rriTy. The notation for the Unes is the same as in Fig. [5] The 
black sohd is m^, the red dashed m^, the blue dotted m^, and 
the green dash-dotted nia- The thick lines correspond to m,, 
and nia at a = 0.05, while the thin lines to a = 0.1. 




FIG. 8: Difference of the screening masses of ao and tt as a 
function of /ifl/m,r and T/rrir: in the extreme case of a = 
0. This quantity indicates the progressive U(1)a restoration 
in the hot and/or dense medium. The masses are denoted 
by m* with asterisk to make clear that these are in-medium 
quantities and different from the vacuum m-^ which is a fixed 
parameter of the model; see Table HI 



tice because the number of lattice sites in the temporal 
direction is severely limited. Instead, one can straight- 
forwardly determine the damping of correlators of the 
fermion-bilinear interpolating fields for the mesons at 
large spatial separations. The exponential decay of the 
correlations is related to the screening mass of the light- 
est mode in the selected channel. In the (P)NJL model, 
this can be found as the pole of the static propagator, 
D{uj = 0,p), in the complex-momentum plane. The re- 
sults are shown in Fig. [T] At T = the real-part and 
screening masses should coincide thanks to the Lorentz 
invariance. In the model calculation, however, this nice 
feature is slightly breached by the three-momentum cut- 
off, but the difference of the two masses in the vacuum 
turns out to be about a few percent at most. Hence the 
apparent cutoff artifacts are reasonably small. 

The masses of the chiral partners become degenerate 
at high temperature, signaling the restoration of chiral 
symmetry. On the other hand, at any (fixed) a 7^ 
the masses of the parity partners, connected by a U(1)a 
rotation, do not converge even at the highest tempera- 
tures considered. From the technical point of view, this 
is a consequence of the simple structure of the inverse 
propagators (|16p : when chiral symmetry is restored, the 
self-energies of the parity partners become equal up to 
a simple rescaling by = (1 — 2a)^. Physically, in re- 
ality, one should expect the coupling a to vary with T 
since it is induced by instantons whose density is expo- 
nentially suppressed at high T j50j. If we consider that 
the U(l)A-breaking interaction strength is proportional 
to the (pure) topological susceptibility, we can infer the 
T-dependence from the lattice data in the pure gauge 
simulation. Instead of doing so, in this work, we pick up 
several values of a. 

A proper way to understand Fig. [7] is thus as follows. 
At T = naturally a is nonzero, and if precise two-color 
simulation data is available for and m^, in princi- 



ple, a can be fixed by the data. We can perform the 
PNJL model calculations using the determined a to go 
to the higher temperature. If we see a reduction of 
and ma toward degenerated and rua, it is a signal 
for the effective U(1)a restoration. We can deduce how 
far a decreases by adjusting a to fit m,, and TOq at each 
temperature. Hence, Fig. [7] is a demonstration for all 
these possible investigations once the two-color simula- 
tion successfully measures the screening masses in good 
precision. 

Finally, in Fig. [SI we plot the difference of screening 
masses of oq and tt as a function of /Lts/TO^ and T/m^. 
Since these two modes do not mix even in presence of 
the diquark condensate, the masses can be calculated 
straightforwardly also in the BEG phase using the in- 
medium propagators (|15p . The choice of the screening 
masses instead of the real-part ones here is motivated 
by the lattice measurement, and also technically favored: 
while the real-part masses are obscured by Landau damp- 
ing in the diquark condensation phase at T 7^ 0, the 
screening masses remain well defined. The results in- 
dicate effective restoration of U(1)a symmetry at high 
temperature and/or chemical potential. While in the T- 
direction the degeneracy more-or-less copies the restora- 
tion of chiral symmetry which is another source of U(1)a 
breaking, in the /^s-direction the convergence of m* and 
TO* is slower as a result of additional U(1)a breaking by 
the diquark condensate. We note that Fig. [8] is the ex- 
treme example of a = 0, which is not likely near the 
phase boundaries, but could be the case in the quark- 
gluon plasma region in view of the lattice data of the 
topological susceptibility In the future, by combin- 
ing the two-color lattice outputs and the PNJL model 
analysis, it would be possible to make a 3D plot of a 
which should approach zero at high T and / or high /i b ■ 
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IV. CONCLUSIONS 

We have adopted the PNJL model as an effective ap- 
proach to two-color QCD and applied it to the case of 
two light quark flavors. This is the simplest case which 
exhibits nontrivial low-energy spectrum due to sponta- 
neous chiral symmetry breaking, and at the same time 
can be simulated by lattice Monte-Carlo techniques. We 
argued that one can fit the parameters of the NJL part 
of the model using physical (three-color) observables, but 
their values have to be rescaled appropriately. Once this 
is done and the quark sector is coupled to the Polyakov 
loop, the model yields locking behavior of chiral and de- 
confinement crossovers as long as /^s is zero. 

We checked older results on the phase structure of the 
two-color NJL model in the plane of T/m^ and /i^/m^, 
and analyzed its modification induced by the coupling 
to the Polyakov loop. The phase transition between the 
normal and superfluid phases is second order for all values 
of the chemical potential considered. In a large part of 
the diquark condensation phase the expectation value of 
the Polyakov loop is small, which resembles quarkyonic 
matter predicted using large- A^c arguments. We carefully 
clarified the realization of quarkyonic matter in the two- 
color system. The baryon number density hb appears 
finite as soon as fis exceeds the mass of the baryonic pion, 
but still the quark contribution to ns is not substantial 
until fj,B surpasses the twice of the dynamical quark mass. 
After then quark degrees of freedom supersede baryons, 
meaning a transition from "superfiuid nuclear matter" 
into "quarkyonic superfluid." 

Our model analysis of the phase diagram is based on 
two simplifying assumptions. The first one is the mean- 
field approximation which treats the system as a gas of 
noninteracting fermionic quasiparticles. This may not 
be quantitatively accurate in some regions of the phase 
diagram such as for fJ-B — "m-n at nonzero temperature 
where the system behaves rather as a dilute Bose gas. 
On the other hand, in the theory of strongly-interacting 
Fermi gases the mean-field approximation is known to 
be reliable at zero temperature. Moreover, the structure 
of the phase diagram concerning diquark condensation 
and chiral symmetry restoration is robust, being a direct 
consequence of the symmetry of two-color QCD. The sec- 
ond assumption is the extrapolation of the gauge part of 
the thermodynamic potential to nonzero baryon chemi- 
cal potential. This has been justified for three-color QCD 
and low baryon chemical potential by a comparison with 
available lattice data and, in fact, is the source of the 
predictive power of the PNJL model. Therefore, we are 
confident about the existence of the quarkyonic super- 
fiuid phase as depicted in Fig . [TJ On the other hand, the 
most recent lattice data (sj suggest that the Polyakov 
loop at a fixed low value of temperature starts to rise 
at /is ^ 3m^, signalling possible deconfinement. This 
certainly presents a challenge to the PNJL model and 
determines the direction of our future research efforts. 

Finally, we studied the dependence of the spectrum 



of collective excitations (scalar and pseudoscalar mesons 
and diquarks) on the strength of the axial anomaly. For 
that sake we introduced a NJL-type interaction with a 
tunable U(l)A-breaking parameter. For all modes we cal- 
culated both the real-part mass and the screening mass, 
which governs the decay of spatial correlators, in order 
to facilitate a direct comparison with lattice simulations. 
Above the chiral restoration/deconfinement temperature 
the masses of the chiral partners become degenerate, as 
expected. At the same time, our results indicate that the 
restoration of U(1)a symmetry in terms of the masses of 
parity partners cannot be hidden by chiral restoration 
unlike the full topological susceptibility, which is good. 
This would naturally incorporate in the model study the 
suppression of instanton effects in a hot and/or dense 
matter. To make this whole argument into a quantita- 
tive predictive framework, further physical input from 
the lattice simulation would be indispensable, e.g. pre- 
cise measurement of ma as a function of T and ^lb- It 
would be also interesting that, since the topological sus- 
ceptibility in the two-color pure gauge theory does not 
cost much, we can compare the inferred a behavior and 
the suppression of the pure topological susceptibility. 
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Appendix: Ginzburg— Landau expansion of the 
thermodynamic potential 

In order to make clear whether the phase transition 
of diquark superfluidity becomes (weakly) first order at 
high /is, we performed the Ginzburg-Landau expansion 
of the thermodynamic potential near the second-order 
transition line and calculated the coefficient of the quartic 
term with respect to A. In the presence of a tricritical 
point and the onset of a first-order phase transition, this 
coefficient would go to zero and then change its sign to 
negative. 

The thermodynamic potential 51 depends only on the 
square of the diquark condensate. In general, 51 depends 
on other condensates, symbolically denoted by x^, as 
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FIG. 9: Ginzburg-Landau quartic coupling along the second- 
order transition line as a function of fiB/m-n. The solid and 
dashed lines are the PNJL and NJL model results, respec- 
tively, with the same input parameters. 



well. In the PNJL model, we have x = W, while in 
the standard NJL model the only other condensate would 
he a. To study the behavior of the diquark condensate 
near the critical temperature, it is most convenient to 
solve the gap equations for the other condensates, i.e., 
drt/dxa = 0. These define Xa implicitly as a function of 
A^, and the thermodynamic potential is then a function 
of A2 solely, that is, fl = rj(A2,Xa(A^))- 

The coefficients of the quadratic and quartic terms in 
the Ginzburg-Landau functional are now determined by 
the first and second total derivatives with respect to A^ , 
evaluated at A = 0. The first derivative vanishes at 
the transition point by means of the gap equation. The 
second derivative defines the effective Ginzburg-Landau 
quartic coupling and is in general expressed as 



A, 



GL 



d(A2)2 



(17) 



9(A2)2 dA^dxa\dxadxbJ dxbdA^' 
The inverse in the second term is assumed in the matrix 



sense. The sign of Aql decides whether the transition is 
of first or second order. In particular in the NJL model, 
this expression acquires a simple form (with an obvious 
notation for the partial derivatives) 



(18) 



In the PNJL model, one has to calculate the 3x3 matrix 
of second partial derivatives. Given the formula for the 
thermodynamic potential ([9]), these are easily evaluated 
explicitly at A = as. 



(27r)3 e^i^ir 

f) r^-n^ f ^^fe Sinh/g^^ 
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f d^] 



Y (cosh/3Ci, + $)2 



(19) 



The NJL limit is recovered by setting $ — > 1— . 

We plot the numerical results in Fig. [9l from which 
we conclude that the Ginzburg-Landau coupling is al- 
ways positive, and that the phase transition hence is al- 
ways second order, within a reasonable chemical poten- 
tial range (below the cutoff A), and in both NJL and 
PNJL models. 
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